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Abstract: We study the factorization of soft and collinear singularities in dimensionally- 
regularized fixed-angle scattering amplitudes in massless gauge theories. Our factorization 
is based on replacing the hard massless partons by light-like Wilson lines, and defining 
gauge-invariant jet and soft functions in dimensional regularization. In this scheme the 
factorized amplitude admits a powerful symmetry: it is invariant under rescaling of indi- 
vidual Wilson-line velocities. This symmetry is broken by cusp singularities in both the 
soft and the eikonal jet functions. We show that the cancellation of these cusp anomalies 
in any multi-leg amplitude imposes all-order constraints on the kinematic dependence of 
the corresponding soft anomalous dimension, relating it to the cusp anomalous dimension. 
For amplitudes with two or three hard partons the solution is unique: the constraints fully 
determine the kinematic dependence of the soft function. For amplitudes with four or more 
hard partons we present a minimal solution where the soft anomalous dimension is a sum 
over colour dipoles, multiplied by the cusp anomalous dimension. In this case additional 
contributions to the soft anomalous dimension at three loops or beyond are not excluded, 
but they are constrained to be functions of conformal cross ratios of kinematic variables. 
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1. Introduction 

Studies of infrared and collinear singularities of fixed-angle scattering amplitudes in mass- 
less gauge theories have a long history (for early results, see for example [1] and [2]), 
and they have led to remarkable insights into the all-order structure of the perturbative 
expansion. 

These studies are not motivated by a purely theoretical interest: in fact, a detailed 
understanding of the long-distance singularity structure of QCD amplitudes is a crucial 
element in predicting high-energy collider cross sections. Indeed, the calculation of observ- 
able cross sections involves intricate cancellations of soft and collinear singularities between 
real and virtual corrections (see for example [3]). Furthermore, with a precise knowledge 
of singularities one can predict dominant higher-order corrections, and in many occasions 
resum certain classes of logarithmically enhanced contributions to all orders [4,5]. 

Our understanding of long-distance singularities is based on the ideas of factorization 
and universality [6] . Fixed-angle scattering amplitudes are functions of Lorentz-invariant 
combinations of external momenta, which are assumed to be uniformly much larger than the 
relevant infrared cutoff, typically given by the scale of confinement; one expects then that 
exchanges of virtual particles with vanishingly small energies, or with vanishing transverse 
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momenta with respect to given external legs, should decouple from hard exchanges, which 
happen at much shorter distances. Such a decoupling is far from apparent in Feynman 
diagram calculations, but it can indeed be proven to all orders in perturbation theory, once 
gauge invariance is enforced by means of appropriate Ward identities [4] . 

The result has a simple structure, explained in detail in Section 2. Briefly, multi- 
leg amplitudes can be organized as vectors, in a vector space spanned by the irreducible 
representations of the gauge group that can be constructed with the given external particles; 
such a vector can be shown to have a factorized structure: each external leg is dressed by 
virtual collinear emissions, building up a colour-singlet 'jet' function; soft gluons exchanged 
at wide angles are assigned to a separate factor, which is a matrix, mixing the available 
colour representations; this matrix of 'soft' functions is then contracted with a vector of 
hard scattering coefficients, which contain no infrared or collinear singularities. 

Given this factorized structure, one may immediately deduce that the various factors 
in the amplitude obey simple evolution equations, embodying the consequences of renor- 
malization group as well as gauge invariance [7] (for reviews of this viewpoint, see [8,9]). 
Evolution equations of this type were derived for the first time for the form factors of 
elementary fields, with a variety of methods [10-12], and were later extended to Wilson 
loops [13-20] and to cross sections and amplitudes of phenomenological interest [21-24]. 
Solving these equations leads to the exponentiation of all infrared and collinear singular- 
ities. The singularities in the exponent are generated by integrals over the scale of the 
running coupling of specific anomalous dimensions, which can be computed order-by-order 
in perturbation theory. A significant step was taken in [25] , where the evolution equation for 
the Sudakov form factor was solved in dimensional regularization. Within this framework, 
infrared and collinear poles are generated by integration over the scale of the L)-dimensional 
version of the running coupling; the results of exponentiation can then be directly com- 
pared with finite-order Feynman diagram calculations; the Landau pole is also regulated 
by dimensional continuation, so that resummed amplitudes can be computed as analytic 
functions of the coupling at a fixed scale and of the dimension of space-time [26]. This 
approach was extended to multi-leg amplitudes in [27], confirming earlier predictions [28] 
on the structure of singularities at NNLO. 

In recent years, the development of novel and advanced techniques for high-order cal- 
culations in QCD and in general gauge theories has stimulated further investigation of 
the exponentiation of infrared and collinear singularities. In particular, great theoreti- 
cal effort has been made to further our understanding of amplitudes in supersymmetric 
gauge theories, and most notably in the maximally supersymmetric J\f = 4 Yang-Mills 
theory. This theory is of special interest for several reasons: it is quantum conformal 
invariant; in the planar limit, it is expected to have a simple, theoretically accessible 
strong-coupling limit, because of its connection with string theory through the AdS-CFT 
correspondence [29]; finally, its amplitudes and anomalous dimensions are of practical rel- 
evance, since they have nontrivial relations with the corresponding quantities in QCD (see 
for example Refs. [30,31], and recent studies of the Regge limit [32,33]). Explicit calcu- 
lations for the four-point function in N = 4 super- Yang-Mills (SYM) theory have led to 
an all-order conjecture [34], suggesting that non-singular terms exponentiate together with 
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infrared and collinear poles, at least for the class of maximally helicity violating (MHV) 
amplitudes. While this conjecture has now been shown to fail starting with the two-loop 
six-point function [35], it is clear that N = 4 SYM perturbative amplitudes must have a 
remarkably simple all-order structure, which may well be brought under full theoretical 
control in the near future. A step in this direction was taken with the discovery of a sur- 
prising duality between scattering amplitudes in momentum space and expectation values 
of Wilson loops taken in an auxiliary coordinate space [36]. Further, remarkable progress 
was made at strong coupling in Ref. [37], where a calculation of the four-point amplitude 
was performed, by adapting string techniques to dimensional regularization. This allowed 
a direct comparison with resummed perturbative calculations, finding an exact matching 
in the structure of long-distance singularities. Recent results in this fast-developing field 
are reviewed in Ref. [38]. 

Most of the calculations just described have been carried out in the planar limit 1 , 
which has special simplifying properties. In this limit, soft contributions to multi-leg 
amplitudes can be further factorized into a product of 'wedges', each one proportional 
to a form factor, since soft exchanges can only take place between adjacent external legs. 
In essence, in the planar limit amplitudes can have only a single colour structure, so that 
the soft anomalous dimension matrix must be proportional to the unit matrix. All soft and 
collinear singularities are then determined by just two colour-diagonal functions: the cusp 
anomalous dimension 7_ftr(a s ) [16-20], and a subleading function G(a s ) [40], responsible for 
single soft or collinear poles. 

It is of great interest to push our understanding of infrared singularities in terms of a 
limited set of anomalous dimensions beyond the planar limit. Indeed, from a theoretical 
point of view, only at non-planar level one begins to see the intricate pattern of colour 
correlations that are characteristic of non-abelian gauge theories: only at this level space- 
time and colour degrees of freedom become explicitly correlated. Furthermore, colour- 
subleading contributions in QCD have important phenomenological effects on resummed 
hadronic cross sections, beginning at the next-to-leading logarithmic order, and the un- 
derstanding of subleading poles would also play an important role in the development of 
infrared and collinear subtraction schemes at higher orders in perturbation theory. Finally, 
recent work [40, 43-48] has highlighted new properties of the functions that generate in- 
frared and collinear enhancements in gauge theory amplitudes and cross sections in the case 
of two hard partons, leading to a better understanding of the process-dependence of soft 
radiation, to the discovery of all-order connections between different physical processes, 
and to the possibility of performing internal resummations of running-coupling corrections 
within the Sudakov exponent. It would be very interesting to extend these studies to 
general colour configurations. 

Soft anomalous dimension matrices for multi-particle scattering have also been inten- 
sively studied in recent years. A complete one- loop calculation for the simplest non-trivial 
case of 2 — ► 2 scattering was carried out originally in [49] . More recently, the calculation was 
reproduced in a physically motivated, dipole-based formalism in [50]: an interesting obser- 

1 An exception is Ref. [55], which studies the leading infrared singularities in subleading colour compo- 
nents of the TV = 4 SYM gluon-gluon scattering amplitude to three-loop order. 
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vation there was that the anomalous dimension matrix for gluon-gluon scattering displays 
an unexplained symmetry relating kinematic invariants with the number of colours N c . A 
different symmetry property was observed by [51], where it was noted that all one-loop 
anomalous dimension matrices are complex symmetric matrices in a suitably chosen or- 
thonormal basis. This property was later explicitly verified with the calculation of the 
matrices for all 2 — > 3 processes at one loop [52,53], and very recently proven [54]. 

Finally, a remarkable result was derived in [56] , where it was shown that soft anomalous 
dimension matrices at two loops, with any number of external legs, are proportional to 
their one-loop value, with the proportionality constant given by the two-loop coefficient 
of the cusp anomalous dimension. This is of course a great reduction in the number 
of possible degrees of freedom, since a priori each matrix element could have acquired 
an independent two-loop correction. The fact that the correlation between colour and 
kinematic dependence in the soft function does not get more complex at two loops as 
compared to one loop, calls for a deeper explanation. At present it is not known whether 
this remarkable property remains valid at higher orders. 

In this paper, we begin to tackle this question. In Section 2 we develop in detail the 
factorization of soft and collinear singularities for fixed-angle scattering amplitudes, fol- 
lowing the approach of Ref. [40]. There are two main differences between our factorization 
and earlier calculations of soft matrices. First, we employ dimensional regularization as 
the unique infrared and collinear regulator: thus, for example, in contrast with Ref. [56] 
we do not tilt the Wilson lines off the light cone to regulate collinear poles. While this 
approach makes explicit loop calculations slightly more delicate, it has the advantage that 
Wilson line correlators are given by pure counterterms to all orders in perturbation theory, 
and they do not depend on any mass scales. Second, instead of using the jet definition 
as the square-root of the Sudakov form factor as in Refs. [27,56], we define each jet Jj by 
introducing a separate auxiliary vector rij, as suggested in early work on Sudakov factoriza- 
tion [57]. This will allow us to conveniently trace the effect of rescaling of the Wilson-line 
velocities. 

Section 3 studies the kinematic dependence of the eikonal functions that enter the 
factorization of multi-parton amplitudes. First, in Section 3.1, we consider eikonal jets, 
and we determine their kinematic dependence to all orders in perturbation theory in terms 
of the cusp anomalous dimension. This simple result follows from the fact that the eikonal 
jet is defined as a correlator of semi-infinite Wilson lines (see (2.5) below) one of which goes 
along the light-like direction defined by the momentum of an external hard parton. Any 
such correlator of semi-infinite Wilson lines is classically invariant under rescaling of any of 
the corresponding velocity vectors (independently of whether they are light-like or not): this 
invariance is a property of the eikonal Feynman rules. In the presence of cusps with light- 
like rays, however, the renormalization procedure breaks this invariance: the counter terms 
include double poles, corresponding to overlapping ultraviolet and collinear singularities, 
along with single poles that carry explicit dependence on the normalization of the light-like 
velocity vectors. Thus, the renormalized correlators, which do retain their invariance under 
rescaling of any non-light-like Wilson-line velocity vector, acquire a dependence on the 
normalization of the light-like ones. The origin of this anomaly is well understood [13-16]: 
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the violation of classical rescaling invariance is governed by the cusp anomalous dimension, 
and we will refer to it as the 'cusp anomaly'. 

In Section 3.2 we extend the analysis to soft gluon functions. To deal with the general 
multi-leg case, we examine combinations of soft and jet correlators where the cusp anomaly 
cancels, so that rescaling invariance must be recovered. We find that this strongly con- 
strains the dependence of the soft anomalous dimension on the kinematics of the scattering 
process, and eventually also on the colour degrees of freedom. 

Section 4 deals with the simple case of amplitudes with only two hard coloured par- 
tons. In Section 4.1 we develop the consequences of the new constraints for the case of the 
Sudakov form factor, and show that the complete dependence of the corresponding eikonal 
function on the kinematics is indeed governed by the cusp anomalous dimension. In Sec- 
tion 4.2 we analyse the kinematic dependence of the partonic jet function and contrast it 
with the eikonal case. 

In Section 5 we return to the case of generic multi-parton fixed-angle scattering am- 
plitudes and study the impact of the new constraints on the soft function. We show that 
while these constraints are insufficient to fully determine the functional dependence on the 
kinematic variables, they admit a remarkably simple solution, where the soft anomalous 
dimension matrix at any order in perturbation theory is proportional to the one- loop result. 
The solution corresponds to a sum over all colour dipoles, which correlate the kinematic 
dependence to the colour degrees of freedom, multiplied by the cusp anomalous dimension. 
This formula is consistent with the result of Ref. [56] at two loops and generalizes it to all 
orders. We also discuss possible sources of further corrections. We conclude in Section 6 
by summarizing our results, while two appendices discuss concrete examples. In Appendix 
A we describe the special case of amplitudes with three hard partons, where the sum-over- 
dipoles formula is the unique solution to the new constraints. Appendix B studies 2^2 
scattering of quarks at one loop, describing the way in which conformal cross ratios are 
formed through a sum over diagrams. 



2. Factorization of fixed-angle scattering amplitudes 

We begin by describing the factorization of a general fixed-angle massless gauge theory 
amplitude into soft, hard and jet functions. We follow the notations of Ref. [40] and 
generalize the definition of the soft function given there to the case of multileg amplitudes. 
The amplitude M. describes the scattering of n hard massless gauge particles (plus any 
number of colour-singlet particles) so it is characterized by n colour indices {ccj}, i = 
1, . . . n, belonging to arbitrary representations of the gauge group. The representation 
content of the amplitude is collectively denoted by [/]. Such a coloured object can be 
decomposed into components by picking a basis of independent colour tensors with the 
same index structure. We denote these tensors by (cl)| q . j, where L = 1, . . . 7V"W and iV"W 
is the number of irreducible representations of the gauge group that can be constructed 
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with the given particles. We write then 



ATM 

M® i} ( Pt /n,a s (n 2 ),e) = ( K / M ,a s ( M 2 ),e) (c L ) K} , (2.1) 

L=l 

with jjL being the renormalization scale and e = 2 — D/2, where D is the dimension of 
space-time. General factorization arguments guarantee that the colour components Ml of 
the amplitude may be written in a factorized form. Following [2,27,40,49,56], we write 

^^a,(,. j ) l£ )=4(ft-M 2 ^) H^(^i, (2p v n ; )2 ,^(/x 2 ) 



?1 

x 



, (2/), • /;,•)- 

J * ( 12~~2 ' 0s ^ )' £ 



ffl ffl 
where the hard function iJ^-, like the amplitude .M^ itself, is a vector in the colour space 

ffl 

described above; the soft function S LK is a matrix in this space, while the jet functions J, 
and Ji do not carry any colour index. A sum over K is assumed on the r.h.s. The soft 
matrix S and the jet functions J and J contain all infrared and collinear singularities of 
the amplitude, while the hard functions Hk are independent of e. Each of the functions 
appearing in Eq. (2.2) is separately gauge invariant and admits an operator definition 
given below. These definitions also clarify the choice of the arguments of each function. 
In particular, we have exhibited here the fact that the eikonal functions S and J depend 
on the dimensionless four-velocities Pi associated with external particles, rather than the 
particle momenta pi. The velocities are defined by scaling the momenta pi according to 
Pi = PiQo/V^, where the magnitude of Qq is unimportant so long as this substitution 2 
is restricted to the eikonal functions. The fixed-angle assumption implies that all scalar 
products Pi • Pj (i / j) are of order 1. 

The definitions of the soft and jet functions all involve Wilson lines, which we write as 



$n(A 2 ,Ai) = Pexp 



r^2 

ig dXn ■ A(Xn) 



(2.3) 



In terms of these Wilson-line operators, one may then define the 'partonic jet' functions 
(for, say, an outgoing quark with momentum p) as 



u{p)j[ (2P 2 n / ,a s (/i 2 ),e) = (p|V(0)$ n (0,-°o)|0>. (2.4) 



The function J represents a transition amplitude connecting the vacuum and a one-particle 
state. The eikonal line <!>„ simulates interactions with fast partons moving in different 



2 As soon as the /?» variables are used in the hard functions H or in the partonic jet functions J, which 
do depend on physical scales, Qo needs to be specified. We shall avoid that, except in Section 4.2 where we 
relate the normalization of the partonic jet to the Sudakov form factor. 
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directions: the direction is arbitrary, but off the light-cone (in order to avoid spurious 
collinear singularities). Since eikonal Feynman rules are invariant under rescalings of the 
eikonal vector n^, and this invariance is not broken by the cusp anomaly for n 2 / 0, J can 
depend on the vectors p and n only through the argument given in Eq. (2.4) 3 . To avoid 
any ambiguity with respect to unitarity phases associated with the first argument of J we 
shall choose n such that p ■ n > 0. Note that this can be done so long as one retains the 
vectors n corresponding to different partons independent of each other. 

The factorization formula (2.2) also requires to introduction of the eikonal approxima- 
tion to the partonic jet J, which we call the 'eikonal jet'. It is defined by 

(2(8 -nS 1 \ 
J ( n2 > a s (ji 2 ),e J = <0| <Moo, 0) <& n (0, -oo) |0> . (2.5) 

Both the partonic jet (2.4) and the eikonal jet (2.5) have infrared divergences, as well 
as collinear divergences associated to their light-like leg; thus, they display double poles 
order-by-order in perturbation theory. The double-pole singularities are however the same, 
since in the infrared region J correctly approximates J: singular contributions to the two 
functions differ only by hard collinear radiation. 

It is important to note that the eikonal jet J depends on the renormalization scale 
only through the coupling: indeed, diagram by diagram J is given by integrals with no 
dimensionful parameter. Such integrals vanish identically in dimensional regularization, 
but since this trivial result involves cancellations between ultraviolet and infrared singular- 
ities, upon renormalization J becomes non-trivial: the contribution of each graph equals 
minus the corresponding ultraviolet counterterm. As a consequence, using a minimal sub- 
traction scheme, the result for J at each order in a s is a sum of poles in e, without 
any non-negative powers. These properties are not special to the jet function, but apply 
to any eikonal function not involving dimensionful parameters, provided it is defined in 
dimensional regularization and in a minimal subtraction scheme. 

The final ingredient in Eq. (2.2) is the soft matrix. It is constructed by taking the 
eikonal approximation for all soft exchanges. Since soft gluons are insensitive to the struc- 
ture of hard collinear emissions, they couple effectively to Wilson lines in the colour rep- 
resentation of the corresponding hard external parton. Such exchanges mix the colour 
components of the amplitude, so one is led to define 

n 

(cl) {q) , } S [ [ ] K ■ /3j,a s (n 2 ), e) = £<°l fi [^(°°. )a fc ^] |0> (c K ) {vk} , (2-6) 

where for simplicity of notation we have defined all eikonal lines as outgoing. Note that in 
our definition we keep all Wilson lines on the light-cone. As a consequence, also the soft 
matrix is a pure counterterm in dimensional regularization and it depends on the renor- 
malization scale only through the coupling; furthermore, the homogeneity of the eikonal 
Feynman rules with respect to rescalings of the eikonal vectors 3\ would suggest that S 



! For later convenience factors of 2 have been introduced into the arguments of the jet functions. 
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can depend on Pi only through homogeneous ratios invariant under such rescalings. As de- 
scribed in [40], this is not true: indeed, rescaling invariance is broken by the cusp anomaly, 
so that the soft matrix acquires nontrivial dependence on the scalar products Pi ■ (3j. This 
observation will be central to our arguments in the rest of the paper. 

The soft matrix, Eq. (2.6), displays both infrared and collinear poles. One must then 
correct the factorization formula in order to avoid double counting of the infrared-collinear 
region for each external leg. This is achieved in Eq. (2.2) by dividing by an eikonal jet Ji 
for each external leg, thus removing from Jj its eikonal part, which is already accounted 
for in S^ K - One may then observe that the ratios Ji/Ji are free of infrared poles, and 
thus contain only single collinear poles at each order in perturbation theory. Similarly, the 
'reduced' soft matrix 

s% faM*) = J^^^ x (2.7) 



[[Ji I — — 2 — ,<w )> e J 

A — 1 N % / 



is free of collinear poles, and contains only infrared singularities originating from soft gluon 
radiation at large angles with respect to all external legs. This means that the effects of 
the cusp anomaly, which is the source of double infrared-collinear poles, must cancel in S. 
More generally, invariance under rescaling of each individual light-like eikonal velocity, 

Pi KiPi , (2.8) 

which is broken separately in S and in must be recovered in their ratio, Eq. (2.7). Indeed, 
in the factorized amplitude (2.2) the dependence on the normalizations of the vectors Pi 
appears only though the eikonal functions contained in Eq. (2.7), so the invariance of the 
amplitude as a whole with respect to such rescaling amounts to invariance of S. The 
immediate consequence is that S can only depend on arguments that are simultaneously 
homogeneus in Pi and in m. Given the different functional dependencies of S and J , this 
can be achieved only if S depends on kinematics only through the variables 4 



Pij ~ m-ni)\P r n 3 f • 



In Section 3.2 we will explore further consequences of this constraint on the functional 
dependence of the reduced soft matrix. 

Finally, it is important to control the ultraviolet behavior of the jet and soft functions 
thus introduced. All these functions are multiplicatively renormalizable [14, 15]; there is 
however an important difference between eikonal operators involving light-like Wilson lines 
and partonic amplitudes. First, as already mentioned, the ultraviolet divergence of eikonal 
operators is directly related to their infrared singularities. Moreover, anomalous dimensions 
of operators involving cusped Wilson lines with light-like segments are themselves divergent, 
due to the overlapping of collinear and ultraviolet poles. These divergences are controlled 



4 Following [28] we keep track of the unitarity phase by writing —f)i ■ /3j = \/3i ■ /3j\ e 1 *^ where Xi 
if i and j are both initial-state partons or are both final-state partons, and X%j = otherwise. 
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by the the cusp anomalous dimension [16-20]. Let us then write down renormalization 
group equations for the various functions defined above. 

The partonic jet J does not involve any light-like Wilson line, and therefore does not 
have a cusp anomaly. Its anomalous dimension is finite, and one may write 



In both cases the index i is kept as a reminder that the jet function J for a given parton 
i carries information not only on the kinematics, but also on the parton spin, flavor and 
colour, while the eikonal jet J depends on the colour representation only. 

For the soft matrix S multiplicative renormalizability must be understood in the sense 
of matrix multiplication [15], and one writes 



V-T^ik (A • Pj, «s( A e) = " I/j (A • Pj, « S G« 2 ), e) Sf K (A • Pj, « s (/" 2 ), e) . (2.12) 



where Tjj will be referred to as the 'soft anomalous dimension'; it is similar to the 'cross 
anomalous dimension' of Ref. [15,22], taken in the limit where all the lines are light-like. 
In the following sections we shall examine the dependence of the anomalous dimensions 
defined above on the kinematic variables as well as on the colour degrees of freedom. 

3. On the kinematic dependence of eikonal functions 

We now discuss the properties of the eikonal functions J and S taking into account their 
gauge invariance, their renormalization group evolution and their independence of any 
dimensionful kinematic scale. By considering the effect of velocity rescaling we deduce that 
the kinematic dependence of these functions is tightly connected with cusp singularities. 
We first illustrate this for the eikonal jet function J, and then we move on to consider the 
central object of our work, the soft anomalous dimension matrix introduced in Eq. (2.6). 

3.1 Explicit solution for the eikonal jet 

Let us consider first the anomalous dimension of the eikonal jet in Eq. (2.11). One observes 
that the homogeneity of eikonal Feynman rules under the rescaling in Eq. (2.8) would forbid 
any dependence on Wi = 2(/3j • rii) 2 /nf, were it not for the cusp singularity. One expects 
then that the full u>i dependence of 7^ should be proportional to the cusp anomalous 
dimension, and this is indeed the case as we now explicitly show. 

Our starting point is Eq. (2.11); in dimensional regularization, the statement that Ji 
is a pure counterterm implies that it can depend on jjl only through the running coupling; 
one may then solve Eq. (2.11) as 




(2.10) 



In contrast, for the eikonal jet Ji we write 




(2.11) 





(3.1) 
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Next, we observe that the eikonal jet must obey an evolution equation of the same form 
as the Sudakov form factor itself, a so-called 'K+G' equation; a similar observation was 
made in Ref. [40] concerning the eikonal approximation to the form factor. Following the 
standard reasoning, one rewrites the anomalous dimension 7^ as a sum of a singular term, 
generated by the cusp singularity, and a residual finite function that contains the kinematic 
dependence. We write then 

Ui {wi,a s (n 2 ,e) ,e) = - ^Gj t (wi,a s ([i 2 ,e)) + 7$ (a s (X 2 , e)) . (3.2) 

Here we have introduced the cusp anomalous dimension j^(a s ), for an eikonal line in the 
representation of parton i, and a remainder function Gj i (wi, a s ). The normalization of 
the singular term on the r.h.s. of Eq. (3.2) is one half of the corresponding term in the 
Sudakov form factor, since the form factor is comprised of two jets 5 . 

Upon inserting Eq. (3.2) into Eq. (3.1), and changing the order of integration, one 
readily arrives at 



Ji (u>;,a s (^ 2 ),e) = exp 



I jff Q^Ko^V)) - \^(*s(e,e)) In |) 



(3.3) 

which is analogous to the expression for the Sudakov form factor, as given in Eq. (4.22) 
below (or in Eq. (2.11) in Ref. [40]), with the physical scale of the form factor, —Q 2 , 
replaced here by the renormalization point fi 2 . The finite function Gj i has no explicit e 
dependence in a minimal subtraction scheme, since Ji is a pure counterterm. 

We are now going to show that Eq. (3.3) can be further simplified, since the dependence 
on the kinematic variable W{ in the function Gj % can be completely solved for. In order to 
do that, we use the results of Ref. [40] for the W{ dependence of the eikonal jet, which is 
given by 

^ln^K*VM = -\jf f^K^sie.e)) . (3.4) 

Clearly, Eq. (3.3) and Eq. (3.4) are compatible only if Gj x (wi,a s ) is a linear function 
of lnwj. Indeed, by taking the derivative of Eq. (3.3) with respect to \nwi, and using 
Eq. (3.4), one gets 

for any fi 2 . Therefore 



d 1 

Wi dw Gjt 0s ^ 2 ' £ ^ = ~ 2 7 * ^ 6 ^ ' 

which we integrate to get 



G Jt (wi, a s ) = -- 7$ (a s ) \n( Wi ) + <5^(a s ) , (3.7) 



5 Indecd an alternative definition of the partonic jet function, which was used for example in Refs. [27,56], 
is based on taking the square root of the Sudakov form factor. 
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where bj i is a constant of integration, free of any kinematic dependence. Using Eq. (3.7) 
in Eq. (3.2) we finally get 

+ \£ 2< fl$(«s(e,e)). (3.8) 
We can now write down our final expression for the eikonal jet, using Eq. (3.3). We obtain 

Ji{Wi,a s {n ),€) =exp< - J _ -^.^(A , e) J - -<y K K ' \a s {\ ,e)J In I 

(3.9) 

where, as anticipated, the kinematic dependence of the eikonal jet is explicitly written, to all 
orders in perturbation theory, in terms of the cusp anomalous dimension. We observe that 
the cusp anomalous dimension simultaneously controls the double poles and the kinematic 
dependence of the single poles. In the following sections we will see that this property 
holds also in the more complex soft functions. 

Returning to the comparison with the Sudakov form factor case (see Eq. (4.22)), we 
now see that the physical scale — Q 2 is replaced here by n 2 Wi = 2/x 2 (/3j • rii) 2 /n 2 . It is 
important to note that Eq. (3.9) can also be expressed as 

Ji(wi,a s (fi 2 ),e) = expj^jf y ^5 Jt (a s (fi 2 6, e)) - ^^\a s (fi 2 9, e)) In (^j-^j j, 

(3.10) 

exhibiting the fact that \x dependence appears only through the D-dimensional running 
coupling. 

Finally we emphasize that the above result for the eikonal jet holds for quarks as 
well as for gluons. In fact, the dependence of Eq. (3.9) on the colour representation of 

(i) 

the parton i appear only though the two functions j K and Sj^. Moreover, the non- 
Abelian exponentiation theorem [41] implies that the colour structure of these functions is 
'maximally non-Abelian'. Up to three loops, this implies, in particular (see e.g. Refs. [18, 
20]) that the cusp anomalous dimension depends on the representation only through an 
overall multiplicative factor, the total colour charge, given by the quadratic Casimir Cj in 
the representation of parton i, 

CJ = ^TS a) Tj a) , (3.11) 

a 

where I is the unit matrix and is a generator in the corresponding representation 6 . 
Casimir scaling, namely the universality of 7^ (a s ) jCi between quarks and gluons, has 
been explicitly verified in recent years by three-loop calculation of the QCD splitting func- 
tions in Ref. [58]. Starting at four loops, however, the colour structure in the exponent may 



6 T^ a ' should be interpreted as follows: for a final-state quark or an initial-state antiquark: t a af) \ for a 
final-state antiquark or an initial-state quark: — tp a ; for a gluon: if ca b- For SU(iV c ) the index a runs from 
1 to N'c — 1, and specifically Cf = Tr(Nc — 1)/N C for quarks and Ca = N c for gluons. In our normalization 
T R = 1/2. 
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not be expressible in terms of quadratic Casimirs. In general, higher Casimir contributions 
do appear in QCD calculations at this order, for example in the QCD beta function [59,60], 
where one finds colour-singlet contributions constructed of traces of products of four gen- 
erators. The potential appearance of such higher Casimir terms in the exponent, despite 
the non-Abelian exponentiation theorem [41], was first observed in Ref. [42], where it was 
proposed to describe the colour structure of the exponent by 'colour connected webs', giv- 
ing a more precise meaning to the notion 'maximally non-Abelian'. Recently it has been 
argued [39], based on different theoretical considerations 7 , that such terms may indeed 
appear in the cusp anomalous dimension starting at four loops. 
Let us therefore write in full generality, 

l^{a s ) = C i ^ K {a s )+^ , (3-12) 

where Cj is given by (3.11), (a s ) = 2a s /ir + . . ., and 7^ = 0(a^.). Note that jk (a s ) 
is a universal function of the coupling, strictly independent of the representation of the 
parton i. This function is known [58] up to three loops in QCD. In contrast, the residual 
term 7^ represents (yet unknown) potential contributions which violate Casimir scaling; it 
depends on the representation in a more complicated way, for example through terms that 
involve irreducible combinations of four colour generators. The particular way in which the 
cusp anomalous dimension depends on the representation will not be important for most 
of what follows, but it will be used in Section 5 for constructing an explicit expression for 
the soft anomalous dimension. 

In a similar way one expects that 8j t of Eq. (3.7) would be proportional to the quadratic 
Casimir at least up to three loops, so we write 

5^( as )^C l 6j(a s )+5 i p, (3.13) 

where Sj(a s ) = a s /ir + 5j? = O(q^), and d is the Casimir operator defined in 
Eq. (3.11). The one-loop result quoted here can be deduced from the calculation in Ref. [40]. 

3.2 Factorization constraints for soft anomalous dimension matrices 

Having established the simple result in Eq. (3.9) for the eikonal jet, where the kinematic 
dependence is determined by the cusp anomalous dimension, one may wonder if the same 
is true for the soft function. In other words, one may ask whether the full dependence of 
Fjj (/3j • f3j, a s (fx 2 ), e) in Eq. (2.12) on • f3j is associated with the cusps. 

One may try to apply the rescaling argument, arguing that if not for the cusp singu- 
larities Vfj should have been invariant with respect to independent rescalings of each 
One realises however that if the number of hard external lines is n > 4, it is possible to 
construct homogeneous conformal cross-ratios such as 

Pm = wmFm (3 - 14) 



7 We thank Juan Maldacena for pointing this out to us. 
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which are inherently invariant with respect to rescalings of each individual velocity, thus 
evading this argument. Kinematic dependence does not necessarily lead to violation of 
the rescaling-invariance property, and therefore might not be associated with the cusp 
singularities. It is important to note, though, that for n = 2, 3 there are no such ratios and 
the argument does hold. One should expect, therefore, that at least for n = 2, 3 the full 
kinematic dependence should be controlled by 7x to all orders. We shall see below that 
this is indeed the case. 

The observation that allows us to make a step forward is that the soft function, for any 
number of legs, can be indirectly constrained by considering the kinematic dependence of 
the reduced soft function, defined in Eq. (2.7). Here the cusp singularity itself cancels out, 
and yet, as we will see, it leaves its trace through the dependence on the kinematics. To 
proceed, consider the renormalization group equation for the reduced soft matrix S, which 
reads 

Pj^^iK {Pij, «s(M 2 ), e) = - rfj (/Hj, a s (p 2 )) S [ j ] K (pij, a s (p 2 ), e) , (3.15) 

where Tfj, in contrast to rfj and jj, is free of singularities. Its invariance with respect 
to scaling of each individual velocity Pi is manifest in its functional dependence on the 
velocities only through the ratios p^, defined in Eq. (2.9). 

Given the definition of the reduced soft matrix in Eq. (2.7), one easily sees that the 
various eikonal anomalous dimensions are related by 

n 

T u(Piji a s) =rfj(A" Pj,a s ,e) - 5/j^7j fc (ti; fe ,a s ,e) , (3.16) 

k=i 

where, as above, wj, = 2(/3/% • rik) 2 /n 2 ,. In words, pole terms must cancel on the right-hand 
side of Eq. (3.16), and the functional dependence on eikonal vectors must be arranged 
so as to reconstruct functions of p^ in order to be consistent with the left-hand side. 
Substituting Eq. (3.8) into Eq. (3.16) we get an explicit expression for T^j. In terms of the 
D-dimensional running coupling, we can write 

n r i 

r /j (pij,a s (p 2 , e)) = rfj ((3i ■ Pj,a s (p 2 , e), e) - 5 U ^ ~ l^J k ( a s(p 2 , 0) 

fc= ! (3-17) 



This immediately implies that 

• off-diagonal terms in T s must be finite, and must depend only on conformal cross 
ratios of the form of pijki in Eq. (3.14), which indeed can readily be turned into ratios 
of pi/s, as defined in Eq. (2.9), for example 

Pijki - ^'fZ k '^\ = (—) V2 e-»^+*«-^> ; (3.18) 

(A ■ wiPj ■ Pi) \PikPjiJ 
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• singular terms in T s must be confined to diagonal matrix elements, and must be 
determined by the cusp anomalous dimension according to 

rf J {f3 i .f3 j ,a s (» 2 ,e),e)=5ijJ2 \ /" ik^ ej) + 0(e°) ; (3.19) 

k=i Jo ^ 

• finite terms in the diagonal matrix elements must include — in addition to terms that 
depend exclusively on conformal cross-ratios, as in Eq. (3.18) — terms with definite 
kinematic dependence on • f3j, which are proportional to 7^, so as to properly 
combine with the ln(iUj) terms in Eq. (3.17). 

To illustrate how these features arise, in Appendix B we perform an explicit one-loop 
calculation of a 2 — > 2 quark-antiquark scattering amplitude. Note in particular that a 
given diagram violates rescaling invariance also in off-diagonal terms, but this violation is 
eliminated upon taking the sum of all diagrams, which is where conformal cross ratios like 
Pijki are formed. This is a consequence of gauge invariance. 

Returning to the general case, in Section 5 we will give an explicit formula that satisfies 
the requirements outlined above. Our goal here is to first formulate the requirements in 
a compact and general way. To this end, let us consider the derivative of Eq. (3.16) (or 
Eq. (3.17)) with respect to ln(iWj). Noting that the Wi dependence appears only through 
the eikonal jet functions, and using the results of Section 3.1, which imply 

d 1 
Wi dm 1Jl a " ^ = 4 7 ^ ^ ' ^ 3 ' 2 °^ 

we obtain a simple result for the Wj-dependence of j, 

Wi — Tjj(p ij ,a s ) = -5 I jw i —jj(w i ,a s ,e) = - -7$ (a s )Su . (3.21) 

This result can be turned into an equation for the dependence of the anomalous dimension 
matrix on its proper arguments, pij, just using the chain rule. Indeed, for any function F 
depending on Wi only through pij, one finds 



We conclude that 



d — 

Yl d\npi T "( pij,a ^ = °' Vi ' / ^ J ' 



5 d^p~ vSlJ {Pt ^ as) = 4 ^ {as) ' yi ' I = J - (3 ' 23) 

As expected, the first equation in (3.23) states that off-diagonal matrix elements of the 
soft anomalous dimension matrix should be logarithmic functions of homogeneous confor- 
mal cross-ratios of pi/s, such as pijki- For diagonal terms, the second equation in (3.23) 
states that inhomogeneous terms are allowed, but they must be proportional to the cusp 
anomalous dimension in the colour representation of parton i. We will explore the conse- 
quences of these constraints in the following sections, beginning with the case of two-parton 
amplitudes. 
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4. Two-parton amplitudes 

In this section we consider in some detail the consequences of the new constraints in the 
simplest case of amplitudes with two hard coloured partons. We choose to analyse in 
Section 4.1 the special case of the spacelike Sudakov form factor of a quark, but the results 
apply, with minor modifications, to any amplitude with two hard partons. In Section 4.2 we 
consider the partonic jet function, which is an important building block in the factorization 
formula, Eq. (2.2), for any amplitude. We use there the results of Section 4.1 to constrain 
the kinematic dependence of the partonic jet, which is significantly more involved than 
that of the eikonal jet considered above. 

4.1 The case of the Sudakov form factor 

Let us consider the implications of the factorization constraints derived above on the sim- 
plest fixed-angle scattering amplitude, the Sudakov form factor. We will see that the 
constraints of Eq. (3.23) lead to a refinement of the results of Ref. [40], since the kinematic 
dependence of the Sudakov soft function can be explicitly determined in terms of the cusp 
anomalous dimension. 

As for any amplitude, our starting point is the factorization formula of Eq. (2.2), which 
here takes the form 



r^Ja^v) = H {2p ^f , a s (/x 2 )) *S(0i- (3 2 , a s (p 2 ), e) 



x 



(2 Pi ■ m) 2 2 
2 J 2~~2 — > a svM ),e 



n hra '-• <«> 



where Q 2 = (p\ + P2) 2 = %Pi ■ P2- For definiteness, we will consider the space-like form 
factor, Q 2 < 0. 

In the case of the form factor, the soft function S is simply the eikonal correlator 
defined by two Wilson lines running along the classical light-like parton trajectories, with 
velocities given by f}\ and fa- Thus 

S (0 1 -p 2 ,a a (jj?),e) = (01^(00,0)^(0,-00)10) . (4.2) 

To determine the kinematic dependence of S, we consider the reduced soft function S, 
which is given by 

5(p 12 ,a^), e ) = nLii7( ^ as(/i2))e) , (4.3) 

where, as before, Wi = 2(/3j • rii) 2 /n 2 and p±2 = {fix • fa) 2 / {W1W2); the latter is specific to 
the spacelike momentum configuration where A12 = so the phase in Eq. (2.9) is absent. 
The reduced soft function S obeys the renormalization group equation 

^ = - 7*(Pi2, ) j , (4.4) 



- 15 - 



which leads to exponentiation. Since S is a pure counterterm, one simply gets 

— f 1 f^ 2 dX 2 / \ 1 

S (p 12 ,a s (p 2 ),e) =exp<-- J 75 (pi 2 , a s (X 2 , e)J ^ , (4.5) 

in analogy with Eq. (3.1) for the eikonal jet. 

Factorization now requires the anomalous dimension 7^ to be a linear function of In p± 2 . 
Indeed, Eq. (3.23) in this case reads 

d7s(li2,«s) 1 / \ ( . 

dlnp 12 = 4 7 ^ } ' (46) 

which integrates to 

ls(pi2,a> s ) = -^1k{oi s ) lnpi2 + %(a s ) , (4.7) 

where 5g(a s ) is introduced as a constant of integration, and does not depend on pi 2 . As 
expected, the dependence of 7^- on the kinematic variable p± 2 is very simple, and is fully 
determined by the cusp anomalous dimension. 

Note that, similarly to what was done for the jet function in Eq. (3.13), we may 
extract from the anomalous dimension 6$ a factor of the Casimir operator of the relevant 
representation, defining 

S^(a a ) = C i 6 s (a a ) + '^\ (4.8) 

where as usual %(a s ) = a s /ir + . . ., 6® = 0(a a ), and 5$ is a universal function of the 
coupling, independent of the colour representation. The one-loop result quoted here will 
be determined in Eq. (4.21) below. 

Using Eq. (4.7) and Eq. (4.5), we can now write down an explicit expression for S, 
where the kinematic dependence is completely solved. We find 

7W , 2^ N f Ill(0l2) ^ d\ 2 ( ,\ 1 f» 2 d\ 2 , ( rx2 , 

S [pi 2 ,a s {n ),€) = exp < — / — lK ^a a {X , e) J - - / _(S^a a (A ,e) 

(4.9) 

Finally, using the definition of S in Eq. (4.3) we obtain an explicit result for the original 
soft function S, 



S • 2 , a s {p ), e) =S (pi2, a s (p ), e) x U i=1 J (wi, a s (p ), e) 
= exp 



where we defined 



$s <m£ ,e) - -j K a s (£ , e) In 



(4.10) 



(<*s) = Sj {a s ) - <% (a a ) , (4.11) 



combining the two constants of integration introduced in Eqs. (3.7) and (4.7). Eq. (4.10) is 
intuitively appealing: the spacelike or timelike nature of the eikonal form factor is associated 
with the explicit logarithm multiplying the cusp anomalous dimension, just as is the case 
for the full form factor. 
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We conclude this section by briefly comparing our results with those reported in 
Ref. [40]. In order to do so, consider the anomalous dimension of the soft function S, 
defined by 

dlnS (ft • ft,a s ,e) 



din /j, 



75 (ft ■ P2,at s ,e) 



(4.12) 



Using Eq. (4.10), we obtain an explicit expression for 75. At the scale p 2 , it reads 
75 (ft ■ P2,a s {fj 2 ,e),e) = - 5 S (a s {p 2 , e)) + \p K (a s (p 2 , e)) In (-ft -ft) 



+ ■ 



1 de 



IK (a s (f,e)) , 



(4.13) 



which is analogous to the result for the anomalous dimension of the eikonal jet in Eq. (3.8). 
One may of course verify the consistency of the various renormalization group equations 
corresponding to Eq. (4.3), observing that 



7s {Pi2,a s ) = 75 (ft • ft,« s ,e) - 7Ji (w 1 ,a s ,e) - -yj 2 (w 2 ,a s ,e) 



(4.14) 



where logarithms of different arguments on the right-hand side nicely combine to form a 
logarithm of the scale-invariant ratio p± 2 , as expected. 

We can compare our final result for 75 in Eq. (4.13) to the result of Ref. [40], where 
the same anomalous dimension is written as 

1 f^ 2 d£ 2 

75 (ft • P2,a s {fJ 2 ,e),e) = - G cik (ft • ft, a s (/i 2 , e)) + - J -^Ik (a s {f, e)) , (4.15) 

defining the eikonal function G c ik. This allows us to relate to 6s, while solving the 
dependence of on the kinematical vectors ft. We find 



G eik (ft • ft,a s 



7: 1k{ois) In (-ft • ft) + 6 s (a s ) . 



(4.16) 



Finally, let us collect the one-loop expressions for the different eikonal functions discussed 
here 8 , showing that the kinematic dependence found in explicit calculations is consistent 
with the general statements we have made. The soft function S was computed at one loop 
in [23] , with the result 



a. 



S(p l -p 2 ,a s ,e) = l-^-C F 



4--m(-ft-ft) 



+ 0(a 2 ). 



(4.17) 



The eikonal jet at this order can be computed by combining Eq. (3.9) with Eq. (3.13), 
obtaining 



Ji (wi,a s ,e) 



4tt 



C f 



1 I 

72 



+ -(l-ln«)) 



+ 0(a 



According to Eq. (4.3), the reduced soft function is then given by 



S (pi2, a s , e) = 1 + ^ C F - 

Z7T e 



1+ 2 ln (^2) 



+ 0(a 2 ) , 



(4.18) 



(4.19) 



8 We choose to work with the quark form factor; in case of the gluon one the overall colour factor Cf 
should simply be replaced by Ca- 



-17- 



which is indeed a function of pi2, consistent with the general expression in Eq. (4.9). Taking 
the logarithmic derivative of Eq. (4.19) one computes the anomalous dimension 

7S(P12, <*s) = ^ C F (l + \ M.m) \ + 0(a 2 s ) . (4.20) 

The p\2 dependence of 7^, as expected, is given by Eq. (4.7), when the one-loop result for 
7^ = 2 Cf a s /iT is used. The leading-order term of the function &g is then given by 

%(a s ) = C F ^ + 0(a 2 ), (4.21) 

while, using Eq. (4.11) and Eq. (3.13), we verify that Ss = 0(a 2 ), in agreement with the 
results of Ref. [40]. 

To summarize, we have shown that the kinematic dependence of all purely eikonal 
functions entering the form factor can be reconstructed using constraints arising from fac- 
torization, together with general properties of pure counterterms in dimensional regulariza- 
tion. The result is that the entire kinematic dependence of these functions is proportional 
to the cusp anomalous dimension. This confirms our initial expectation, based on the fact 
that violation of rescaling invariance in such functions can only be introduced by the cusp 
anomaly. In Section 5 we will return to the generic case of multi-leg fixed-angle scatter- 
ing amplitudes, and work out the consequences of our constraints. Before that, however, 
let us briefly contrast our findings concerning eikonal functions with the case of partonic 
amplitudes. Again, we choose the simplest example, that of the partonic jet. 

4.2 Constraining the partonic jet function 

We have seen that the kinematic dependence of eikonal functions is remarkably simple, to 
all-order in perturbation theory. The singularity structure of a partonic amplitude resem- 
bles that of the corresponding eikonal function. Nevertheless, because partonic amplitudes 
do depend on dimensionful kinematic variables and have infrared singularities as well as 
distinct ultraviolet renormalization properties, they do not admit similarly simple all-order 
expressions. 

An exception is the Sudakov form factor T(Q 2 , e) itself, which does not get renormalized 
(dT(Q 2 , e)/ln/x = 0), so that its entire perturbative expansion in dimensional regularization 
is determined by its infrared singularities [25,40], according to 

T(Q 2 , e) = exp I \ f~ Ql ^ [ G (-1, a s (f, e), e) - \ lK (a^ 2 , e)) In (^) ] J . 

(4.22) 

Indeed, in this case the eikonal function in Eq. (4.10) is very similar to the full form factor, 
and apart from the different scales, the only qualitative difference is the appearance of 
non- negative powers of e in the function G (— 1, a s , e). 

In a general scattering amplitude, which does get renormalized, the dependence on the 
renormalization point fi 2 and on the kinematic variables is of course distinct, and neither of 
them is associated exclusively with the cusp anomalous dimension. The simplest example 
is that of the partonic jet, which we now consider. 
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Using the factorization formula for the Sudakov form factor, Eq. (4.1), together with 
Eqs. (4.22), (3.9) and (4.10), one can directly determine the product 



(2gj • rij) 2 

-"I 9 1 2 9 ' 1 

V A* < A* 2 
where we defined 



J 



(2pi • «i) 5 



with 



cxp 



-,a s ,e 



+ 



1 r° a dA 2 



4 v a 2 

1 ^ 2 dA 2 



4 7o A 2 L 2 



ip 

i=l 



(2p» • n^) 2 

2 9 ) 5 ^ 



IF 

i=l 



(2pi ■ m 



9 9 ) )f I ) 



-^MA 2 ,e)) Inf^) + G(-l,a s (A 2 ,e),e) 



A 2 , 



n 2 A 2 



(4.23) 



(4.24) 



5j(a s , e) = %(a s ) + G(-l, a s , e) 



(4.25) 



where is the function introduced in Eq. (4.7). Note that 5j(a s , e), in contrast to Sj(a s ) 
in Eq. (3.7), depends on e explicitly, and it does have non-negative powers of e coming from 
G(-l,a s ,e). 

We now observe that in Eq. (4.24) all infrared singularities emerge from the A 2 —> 
limit of the integral over the running coupling in the last line. These singularities must 
all reside in the partonic jet function J, while the renormalized hard coefficient function 
H is finite. This observation, together with the renormalization group equation for J, 
Eq. (2.10), and the fact that J does not depend explicitly on Q 2 , implies that J takes the 
form 



J(4,a s (/i 2 ),e 



exp i hj (a s (p 2 n )) + \ j 2 ^"7J («s(A 2 )) 



(4.26) 



+ 



d\ 2 



-7^(a s (A 2 ,e))ln(g]+i<5 J (a s (A 2 ,e),e) 



where we defined p 2 = (2p • n) 2 /n 2 for brevity. All the functions appearing here are finite 
as e — > 0, so also here singularities are generated only through the integration over the 
running coupling in the last line. The function 5j(a s (X 2 ), e) depends on e explicitly, while 
7k, 7j and hj do not. As far as hj is concerned, this last statement is not obvious a priori, 
and it will be proven to all orders below. 

The structure of the result in Eq. (4.26) for the partonic jet is intuitively clear: infrared 
and collinear singularities are generated by the integration over the scale of the running 
coupling in the second line. These terms in the exponent are similar to the expression found 
for the eikonal jet, Eq. (3.9), and indeed the double poles of the two expressions match, 
as they must. Single poles, on the other hand, are different, on account of hard collinear 
radiation, which is not correctly approximated by J . Further dependence on the hard 
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scale p\ arises in the first line of Eq. (4.26), due to the non-trivial ultraviolet behaviour 
of J, which is dictated by 77. It is natural in this case to start the renormalization group 
evolution at = p 2 , where the function hj gives the initial condition. 

Once again, we may compare our all-order expression with one-loop results. One may 
start with the well-known result [25] for the function G, 



G 



Q 2 



2 ; ^-S; ^ 



7T 



C F 



1 3 

_ + __ e 



7T 

12 



-4 + 0(e 2 ) - 



1 



+ <D{al). 



Using Eq. (4.25) and Eq. (4.7) one then finds 



7T 



5j(a s ,e) = — C F - - e — -4 + 0(e 2 ) + 0(a 



12 



The renormalization of the jet, on the other hand, yields 

Ij {a s ) = -\C F '^ + 0(al). 

4 7T 



(4.27) 
(4.28) 

(4.29) 



Inserting the one-loop results of Eqs. (3.12), (4.28) and (4.29) into our general expression, 
Eq. (4.26), and expanding to 0(a s ), we obtain finally 



J 



47T 



In 



--In 2 



+ ln 



7T 



- i --- + 0( e )[> +0(« 



(4.30) 



Comparing this result to the explicit one-loop calculation, reported in Section 3 of Ref. [40], 
we find full consistency and determine the finite coefficient hj, 



a 



7T 



7T 



hj(a,) = —lC F - + — + (D(aj) 



12 



(4.31) 



We conclude this section by comparing our parametrization of the partonic jet, Eq. (4.26), 
with the results of Ref. [40]. This will then be used to relate the function hj(a s ), which 
sets the normalization of the jet function to the hard function appearing in the factorized 
the Sudakov form factor. We begin by considering the differential equation that controls 
the dependence of J on p 2 . This is, once again, a l K + G" equation, and can be written 
as [40] 



91nJ(p 2 /p 2 ,a s (/i 2 ),e) 1 ( p\ _ 
^2 = 2 Q \ )' e 



l£ 2 ^K(a s (X 2 ,e)) , (4.32) 



where the cusp-related singularity was explicitly extracted. The renormalization group 
equation for J, Eq. (2.10), implies that 



d 51n J (p 2 /// 2 , a s , e) d din J (p 2 /^ 2 , a s , e) 



din fi 2 



dlnpl 



d\npl 



<iln [i 2 



0. 



(4.33) 
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Applying this to Eq. (4.32) in turn implies that the scale dependence of Q is controlled by 
the cusp anomalous dimension, according to 

"' g(^,a s (» 2 ),e) =\ lK {a s ), (4.34) 



cHn/i 2 \p, 
or, upon integration, 

^,a s (^),ej =g(l,a s (p 2 n ),e)- l j 2 — lK {a s (\\e)) . (4.35) 

Now, using our general expression for the partonic jet in Eq. (4.26), we can readily 
interpret the function Q in (4.32) in terms of the anomalous dimensions 7j and 5j, and of 
the finite function h j. We find 

Q (1, a s , e) = f3(a s , e) -^—hj(a s ) + jj(a s ) + ^<5j(a s , e) , (4.36) 
aa s 2 

where the first term is proportional to the (3 function, (3(a s , e) = da s {^L 2 , ej/dln/j,, so that 
it starts contributing at 0(a 2 ). 

Using the definition of 5j(a s , e), given in Eq. (4.25), and setting 9 fi 2 = p^, we may now 
obtain an expression for the function G appearing in the Sudakov form factor, Eq. (4.22), 
which we can compare to the results of Ref . [40] . We find 

G (-1, a s , e) = -2 (3(a s , e)^- hj(a s ) - %(a s ) - 2 7j (a s ) + 2 0(1, a s ,e). (4.37) 



Eq. (4.37) can be compared with the result of Section 4 in Ref. [40], which reads 

d 2 
G(-l,a s ,e) = P(a s ,e) — lnH(-l,Wi,a s ) - ~fs(Pi2,a s ) - 2^j{a s ) + s ^Q i (w u a Sl e) , 



(4.38) 

where we have chosen a specific common normalization for the velocities, pi = (3iQo/\^2, 
with Qo = \/—Q 2 , corresponding to (2pi • rii) 2 /n 2 = —Q 2 Wi, and (3\ • 02 = — 1 5 so that 
p±2 = l/(wi ^2)- Using now our result for jg(pi2, ct s ), Eq. (4.7), we see that Eq. (4.38) 
takes the form 

d 1 

G(-l,a s ,e) = (3(a s ,e)- — lnH(-l,Wi,a s ) + 7 7x(a s ) hi(iui w 2 ) 

da s 4 

2 (4.39) 
- %(« s ) - 27j(a s ) + S ^ j Gi{w i ,a s ,e) . 

i=i 

Eq. (4.39) holds for arbitrary w\ and u>2, implying that the dependence of the r.h.s. on 
these parameters through the explicit term, and through the functions H and Q, must 
cancel out. 



9 Note that in setting (2p • n) 2 /n 2 = /i 2 we assume that the auxiliary vector n is timelike, n 2 > 0. A 
spacelike n is also possible of course, and this choice would be reflected in replacing Q(l,a s ,e) below by 
g(-l,a s ,e). 
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Comparing Eqs. (4.39) and (4.37) we obtain 

-2/3(a s ,e) -^—hj(a s ) + 2Q(l,a s ,e) = \~ jK (a s ) \n(wiw 2 ) 
aa s 4 

d 2 (4.40) 

1=1 

We already know that the dependence of the r.h.s. on the Wi cancels out, so we are allowed 
to use this equality for any u>i. Picking w± = u>2 = 1, we obtain 

lij(a,) = ~lDff(-l ) in = l,a 1 ), (4.41) 

which determines the function hj in terms of the hard function H of the Sudakov form 
factor. Eq. (4.41) implies in particular that hj(a s ) does not carry any explicit dependence 
on e, but depends only on the coupling, as anticipated. 

5. Multi-parton amplitudes: the soft anomalous dimension matrix 

Let us now consider the general case of multi-parton amplitudes and examine the conse- 
quences of our new constraints on the all-order structure of the soft anomalous dimension 
matrix. We first note that the constraint of Eq. (3.23), which can be written as 

g ln( \ T mn {pij,a s ) = - Ik ( a s) $mn , Vi , (5.1) 

is an equality between colour matrices, valid in any basis. We will then proceed to work 
with the colour generator notation, as in [28], without specifying a basis; consequently, we 
will drop the explicit matrix indices M and N in the following. 

Next, we observe that Eq. (5.1) effectively relates the colour structure of the soft 

matrix, which is a priori very complicated, to the much simpler colour structure on the 

(i) 

right hand side. Clearly the particular way in which the cusp anomalous dimension 7^- 
depends on the colour representation of parton i becomes important. Following Eq. (3.12) 
we can write 



Ed _,-o . ,1 
.^.dH^) T = 4 

3,3^ J 



Vi, (5.2) 



where in the first term the dependence on the representation is explicit, while in the second 
it is implicit. Using the linearity of these equations we can obviously write the general 
solution as a superposition of two functions 

r^ = r| c . + r| a (5.3) 

which are, respectively, solutions of the equations 

£ r lc (PV,«.) = \ (£ T« T[A % (a.) , , (5.4) 



Yl d\n(pi ) r ^°- ° s) = 4 ^ ( " s) ' ^ ' (5 ' 5) 
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Here Q.C. and H.C. stand for Quadratic Casimir and Higher (order) Casimir, respectively, 
reflecting the different group theoretical structure of the two contributions . In the first 
equation we exhibited the dependence on the colour representation of parton i, using 
Eq. (3.11). In the following we will focus on determining Tq c leaving aside r§ c = C(af ), 
which is driven by yet unknown corrections to the cusp anomalous dimension, that are not 
proportional to C{. For simplicity of the notation we henceforth drop the subscript Q.C. 
A solution for V s , obeying Eq. (5.4), is given by 



— 1 n - 

r 5 ( Pij , a.) = -i jk (a s ) £ £ H Pij ) £ if if 

ansatz o — * — *■ — * J 

* =1 i>jV« a 
1 n 



(5.6) 



=1 a 



where and were defined in Eq. (3.12) and in Eq. (4.8), respectively. Note that in 
Eq. (5.6) the first term, which couples each pair of partons into a colour dipole, carries the 
entire dependence on kinematics, correlating it with the colour structure, while the second 
term is independent of kinematics and is proportional to the unit matrix in colour space. 
Such a term, in a different factorization scheme, such as the one adopted in Refs. [27,56], 
could be associated with jet functions. In our case, since we start with specific operator 
definitions for the jets, we may find leftover colour-diagonal contributions in r . Note also 
that Eq. (5.6), which is valid for a general n-parton amplitude, reduces to Eq. (4.7) in the 
n = 2 case of the Sudakov form factor. This relation was used in determining the second 
term in Eq. (5.6), which is obviously not constrained by Eq. (5.4). 

The explicit calculation of Ref. [56] at two-loops has established that Eq. (5.6) is the 
full answer to this order. In particular, the soft anomalous dimension does not contain 
correlations involving generators of three different partons, such as f a b c if T^Tj^ (see 
Eq. (5.13) below), despite the fact that single poles in e at two loops are known to contain 
such terms. Ref. [56] has demonstrated that these terms are generated at two loops only 
upon expansion of the product of the soft and hard functions. 

To verify that Eq. (5.6) satisfies Eq. (5.4), let us take a derivative with respect to 
ln(pjj), for specific partons i and j, with i / j. We find 

where we used the fact that Yl a if T^ = S a T T^ . Next we sum over j for fixed i, 
as in the l.h.s. of (5.4), and we get 

V {pt ^ as) - -Uk (a ) V VT (fl) T W 

J 3,&i a (5.8) 

= -^(«,)E T , (a) H")' 
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which coincides with the r.h.s. of Eq. (5.4), as required. The second line of Eq. (5.8) follows 
from colour conservation, which is expressed in the colour generator notation simply by 



n 

E T * (a) = - ( 5 - 9 ) 

Clearly, our ansatz does not in general provide a unique solution. Indeed, Eq. (5.4) is a 
set of n linear differential equations in the variables ln(pjj), while the number of variables 
is quadratic in n, so there may be contributions beyond Eq. (5.6). To clarify the nature of 
possible corrections to our ansatz, let us define 

+ A s MN ( Pij ,a s ) (5.10) 

ansatz 

and summarize our knowledge of A^^. First, as discussed in Section 3.2, and shown 
explicitly in Section 4.1 and Appendix A, A^^y vanishes identically for n = 2,3. Thus it 
may contribute only to n > 4 parton amplitudes, starting at three loops. Second, according 
to Eq. (5.4), it should satisfy the homogeneous equation: 

S d^) A ^ N (/% '' as) = ° V *' M ' N ■ (5 - 11} 

Eq. (5.11) is solved by any function that depends on pij only through conformal cross 
ratios such as Pijki, defined in Eq. (3. 14) 10 . As an example, for four external partons there 
are only two independent conformal cross ratios (note that I/V1423 = P1234 P1342); and the 
general solution can be written as 

Af IN = Aff N (P1234, P1342, a s ) ■ (5-12) 

Note that Eq. (5.11) does not restrict the functional dependence of Af^jy on conformal 
cross ratios, allowing in particular non-linear dependence. Interesting examples, still for 
four partons, are 

^[f? = E E 1 T " T * T < ln (Pi**) ln (PMj) ln (Pujk) , (5.13) 

j,k,l a,b,c 

and 

S m = E E d ^ T J T k T i ln ' (PW) 1^ (Piklj) In 2 (pujk) , (5.14) 

j,k,l a,b,c 

where in both equations the sum over partons is understood to exclude identical indices. 
As already mentioned, is known to appear 11 as part of the 0(1/ e) coefficient in two 

10 As explained in Section 3.2 and illustrated in Appendix B, conformal cross ratios appear in V s also in 
(5.6), starting at one- loop. 

xl In Ref. [56] Hj 2 ^ is written in terms of the three kinematic invariants of a four-leg amplitude, assuming 
momentum conservation. It is straightforward to show that the same expression can be written in terms of 
conformal cross ratios, where momentum conservation is not enforced. 



^MN {piji a s) — F MN (pij,Ot s ) 
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loops amplitudes. Nevertheless, as was shown in Ref. [56], it does not appear in the soft 
anomalous dimension at this order. 

For multi-leg amplitudes with more than four legs, the space of possible solutions to 
the available constraints increases further, as there is an increasing number of conformal 
cross ratios. It remains for future work to decide whether contributions beyond the ansatz 
of Eq. (5.6) do indeed appear. 

We conclude by working out the consequences of Eq. (5.6) for the anomalous dimension 
of the original soft matrix S. Using Eq. (3.16), we write 

_ n 

T/J (ft • Pj, «s(M 2 ), e ) I ansa t z = T U (Ptf> a s{^)) + hj ^ Ui (^;, « s (^ 2 ), e) , 

i=l 

(5.15) 

where, as before, to* = 2 (ft • rij) 2 /n 2 . Substituting here Eq. (5.6), together with the 
expression in Eq. (3.8) for the eikonal jet anomalous dimension, we obtain, after some 
algebra 

n 



V s (ft • ft, a s (^), e) | ansatz = -\ % (a s (^)) ^ ^ In (ft • ft e 1 ^) £ T<°>T 



(o) 
3 



1 jzj^i 

2 



+ (/r))+T / — ^A-(n,,(A-.r)) J^T^T? , 

j=l a 

(5.16) 

where, as in Eq. (4.11), 

S s (a s ) = Sj(a s ) - %(a a ). (5.17) 

The manipulation we applied in order to write Eq. (5.15) in its final form is similar to the 
one we used in Eq. (5.8), but in the reverse order: here one first rewrites the colour-diagonal 
terms that depend on ln(wj), which originate in the eikonal jet contributions, in terms of 
a sum over all other partons j, using = — Ylj^i^j^'i tnen one combines these terms 
with the ln(pij) terms in T s ; finally, using Eq. (2.9), one observes that V s depends only on 
ft • ft, as it must. It is straightforward to check that Eq. (5.16) reduces to Eq. (4.13) for 
the Sudakov form factor. Once again, the second line in Eq. (5.16) is proportional to the 
identity matrix in colour space, and could be associated with jets in a different factorization 
scheme. Indeed, the square bracket gives just one half of the kinematics-independent terms 
in the soft anomalous dimension 75 for the Sudakov form factor, Eq. (4.13). This is the 
contribution that needs to be combined with our jets in order to reconstruct the choice of 
Refs. [27,56], where jets are defined as square roots of the form factor. 

Eq. (5.16) can be compared with the direct computation of the soft anomalous dimen- 
sion at one loop. For the latter we express the sum over all diagrams as 

n 

S (ft ■ ft, a„ e) = 1 + g £ $ 0* • & > £ ) E + °^s) ■ ( 5 - 18 ) 

where the basic one-loop integral, stripped of colour factors, is 

l\f (ft • ft , e ) = I - I ln(ft • ft e i7rA -) . (5.19) 
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The anomalous dimension corresponding to Eq. (5.18) is 

T s (Pi ■ ft, a., e) = E ^ " HPi ' 0j ^ £ + ° (a ' } ' (5 - 20) 

which coincides with Eq. (5.16) upon substituting there the one-loop values of the anoma- 
lous dimensions, and upon using again colour conservation 12 for the l/e term. 

6. Conclusions 

We have explored the singularity structure of fixed-angle scattering amplitudes in massless 
gauge theories using their factorization into gauge-invariant hard, soft and jet functions. 
The factorization formula we use, Eq. (2.2), is based exclusively on dimensional regulari- 
sation. It has the advantage that Wilson line correlators, such as the soft function S, are 
pure counterterms to all orders in perturbation theory, and they do not depend on any 
mass scales. 

Our central observation is the presence of a symmetry under rescaling of the velocities 
that define the eikonal functions S and J, Eq. (2.8); this symmetry is built into the eikonal 
Feynman rules, but it is broken for any eikonal function involving light-like segments, due 
to the cusp anomaly. Multi-parton scattering amplitudes depend on a specific combination 
of eikonal soft and jet functions, the reduced soft function of Eq. (2.7), where the cusp 
anomaly cancels and the rescaling symmetry is restored. We have shown that the specific 
way in which this symmetry is broken and then restored imposes tight constraints on the 
functional dependence of these functions on the kinematic variables, as well as on the 
colour variables. The all-order constraints on the reduced soft function are summarized by 
Eq. (3.23). 

For purely eikonal functions which do not depend on any kinematic scale and are 
affected by the cusp anomaly, we observe that the complete kinematic dependence of the 
single pole terms in the exponent is tightly connected to the double poles, and they are 
both associated with the cusp anomaly. This is easy to see for eikonal functions involving 
two rays, such as the eikonal jet or the eikonal version of the Sudakov form factor, but the 
result turns out to be more general, and through Eq. (3.23) it carries over to the multi- 
leg case. This is contrasted with partonic amplitudes, such as the partonic jet function, 
Eq. (4.26), that have more complicated dependence on the kinematic variables, since they 
depend on dimensionful scales and have non-trivial renormalization properties. 

Our conclusions concerning large-angle soft singularities based on Eq. (3.23) can be 
summarized as follows: for amplitudes involving two or three hard external partons (and 



2 Note that in explicit calculations in a given colour basis the fact that 



y^ T (a) T (q) 

3, i¥=i 



= -as 



iOMN , 



which is a statement of gauge invariance, provides a useful consistency check on the calculation of the colour 
factors. A simple example is provided in Eq. (B.2): one can verify that the sum of the three matrices there 
yields —Cf times the unit matrix, as it should. 



- 26 - 



any number of non-coloured particles) the kinematic dependence of the soft function is 
completely determined, and it is governed, to all orders in perturbation theory, by the cusp 
anomalous dimension. We emphasize that this conclusion does not depend on the way in 
which the cusp anomalous dimension itself depends on the colour representation of a given 
parton. The results for the anomalous dimensions of the reduced soft matrices in these 
cases are given in Eq. (4.7) and Eq. (A. 3). 

For amplitudes involving four or more hard external partons the available constraints 
on the soft function are insufficient to uniquely determine the kinematic dependence. They 
nevertheless relate it to the cusp anomalous dimension. Considering the component of the 
cusp anomalous dimension that is proportional to the quadratic Casimir, we have found 
a minimal solution, Eq. (5.6), for the reduced soft anomalous dimension or, equivalently, 
Eq. (5.16) for the original soft anomalous dimension. This solution, which is valid for 
any number of legs, satisfies the all-order constraints and is consistent with all explicit 
calculations available to date. Eq. (5.6) is written as a sum over colour dipoles: it correlates 
the kinematic dependence with the dependence on the colour variables for each pair of 
hard partons in the very same way as one-loop diagrams do. In contrast with the two- 
and three-leg cases, for four or more legs the solution is not unique: the soft anomalous 
dimension may receive additional contributions (at three loops or beyond) which however 
depend on kinematic variables exclusively through conformal cross ratios. We point out 
that the absence of such contributions at two loops is non trivial: terms that satisfy this 
requirement, displayed in Eq. (5.13), do appear in two-loop diagrams at C(l/e), but the 
calculation of Ref. [56] proves that they do not contribute to the soft anomalous dimension 
at this order. This suggests that Eq. (5.6) is in fact the full answer to any loop order, aside 
from corrections (which in turn satisfy Eq. (5.5)) that are induced by possible contributions 
of higher-order Casimir operators to the cusp anomalous dimension itself, which may appear 
at four loops or beyond. 

To summarize, we took here a step towards the understanding of the infrared singular- 
ity structure of gauge theory amplitudes to all orders in perturbation theory. We did so by 
identifying a rescaling symmetry of eikonal functions appearing in the factorization of the 
amplitudes, which is broken by the cusp anomaly, and then restored in a specific way. Using 
this observation, we derived all-order constraints on soft anomalous dimension matrices, 
that are valid for any number of external partons. We then studied the consequences of 
these constraints and established the complete solution for the soft anomalous dimension 
for amplitudes involving two or three partons, and a minimal solution for four partons or 
more. By formulating our results for the jet and soft functions in terms of a few universal 
anomalous dimensions, which depend only on the dimensionally-regularized coupling, and 
not on kinematics, or on a colour basis, we significantly increase the predictive power of the 
factorization formula, providing powerful checks on multi-loop calculations and a better 
starting point for soft gluon resummation. 
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Note added 

A few hours before the submission of our paper to the arXiv, T. Becher and M. Neubert 
published an independent study [61], proposing an ansatz for the soft anomalous dimension 
matrix which is essentially equivalent to our Eq. (5.16). 
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A. The case of amplitudes with three partons 



It is well known that the colour structure of the soft function for amplitudes with three 
hard partons is trivial: the soft matrix is proportional to the identity matrix. Here we 
briefly explain this property and observe another special property of the soft function in 
this case: similarly to the n = 2 case, discussed in Section 4.1, for n = 3 the contraints of 
Eq. (3.23) completely determine the dependence on the kinematics. 

The reduced soft matrix, defined in Eq. (2.7), in this case is given by 



C[/] / N 

•^MN \Pl2, P23, P31, Ot s , e) 



S MN (A " #2, 02 ■ 03, 03-Pi, a s , e) 



(A.l) 

We now note that Eq. (3.23) yields three independent differential equations, which are 
linear in ln(pjj). Explicitly, they are 







<91n(pi2) 
d 



+ 
+ 







<9hi(/)3i) 
d 



rffjv (pi2, P23, P3l,a s ) = jSmnIk* ( a s) » 



<91n(p 2 3) dln(p l2 ) 

d d 
+ 



<91n(p 



23 J 



<91n(p 3 i) 



T MN (P12, P23, P31,a s ) 
rffjv (pi2, P23, P31,a S ) 



■^SmniP (a 8 ) , 



(A.2) 



1 



Having three independent linear equations in the three variables there is a unique solution. 
It is given by 



rffAT (pij, a s) 



4 ] - 7? - 7?) ln(p 12 ) + ( 7 ? - 7g } - 7^) ln(p 23 ) 



+ 



(7f-7S?»-7?)MPi3) 



+ 



4 1} + 4 2) + 4 3) 



(A.3) 



<5mjv , 
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where a constant term was added, as in Eq. (5.6). Note that we have shown that 
is proportional to the unit matrix in colour space without specifying the representation 
of the partons. This result is completely general, and it generalizes the previously known 
two-loop result to all orders in perturbation theory. Just as in the case of n = 2, the 
entire kinematic dependence is controlled by the cusp anomalous dimension. As explained 
in Sections 3.2 and 5, the fact that uniqueness can be established for n = 2, 3 but not for 
n > 4 is related to the absence of conformal cross ratios of the form of Eq. (3.18). 

In Eq. (A. 3), the dependence on the representation of the various partons is implicit, 
appearing through the functions 7^ (a s ) and 5^\a s ). As discussed in Section 3.1, these 
may include higher-order corrections that are not proportional to the quadratic Casimir. In 
this respect the result of Eq. (A. 3) goes beyond the ansatz of Eq. (5.6), which considers only 
terms that are associated with the quadratic Casimir contributions to 7^ in Eq. (3.12). 
Of course, ignoring 7^ in Eq. (3.12) one recovers Eq. (5.6). Upon substituting Eq. (5.6) 
into Eq. (A. 2), these three equations yield 

£ T (o) £ T (a) +T («)j =_C U 
a 

^tW(tW + T w)=-c 2i (a . 4) 

a 
a 

which are satisfied owing to colour conservation, Eq. (5.9). Eq. (A. 4) also implies that 



(A.5) 



which is consistent with the observation that all colour factors entering Eq. (5.6) in this 
case are proportional to the unit matrix. Finally, the explicit sum-over-dipoles solution to 
Eq. (A. 2) takes the form 



(a) T (o) 
1 1 2 


= c 3 


-c x 


-c 2 , 


O)rp(a) 
2 L 3 




-c 2 


-c 3 , 


(a) T (a) 
3 x l 


= c 2 


-c x 


-c 3 , 



r Q.C. i/Hj,OCs) = -0% (<*a) 



(C 3 - Ci - C 2 ) Hp 12 ) + (Ci - C 2 - C 3 ) ln(p 23 ) 



+ (C 2 -C 3 -C 1 )\n(p 13 ) 



(A.6) 

+ (C\ + C 2 + C 3 ) . 



B. The case of qq — > qq scattering at one loop 

The four-parton amplitude qq — > qq provides a simple example where the colour matrix 
structure is non-trivial. We perform the calculation along the lines of Sec. IV of Ref. [56], 
but factorize the amplitude as in Eq. (2.2), using light-like Wilson lines. 
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The velocities are defined by: 



g(/3i) + 9(ft) ->?(/%) + ?(&). 

It is convenient to set 0i = — v±, (3 2 = —V2-,03 = v 3 and = V4, so that all the scalar 
products Vi-Vj > 0. Note that below we will formally treat the four velocities as independent 
variables, choosing not to enforce explicitly momentum conservation. 
Following Ref. [56], we pick the colour basis 



C2 = #13 #24 



(B.l) 



and we use the convention of Eq. (2.6) to write the result in a matrix form. There are six 
one-loop diagrams altogether, and for each one of them the loop integral yields Eq. (5.19). 
Computing the colour factors in the chosen basis and summing up the contributions of the 
six diagrams according to Eq. (5.18) we get 



S (vi ■ Vj,a s ,e 



+ 



+ 




2 1 1 

— ln(vi -v 3 ) ln(t> 2 • V4) 



1 



ln(?;i ■ V4) ln(^2 • V3) 



-5- (ln(vi • v 2 ) + ivr) (ln(w 3 • V4) + iir 



(B.2) 



+ 0(a 2 s ). 



The corresponding one-loop soft anomalous dimension matrix is then 



V s (vi -Vj,a s ,e) = —C F 

7T 



-\ + In ((vi • v 2 ) (v 3 ■ V4)) 

-\ +ln((vi -v 3 )(v 2 -V4)) 



+ 



a s 
2vr 



— In 

N c 



(vi-V4) (V2-V3) 
(v 1 -v 3 ) (v 2 -v 4 ) 



In 



(vi-V2) (vtyvi) 

(V!-V4) (v 2 -V 3 ) 



- 2vri 



TT 




Finite terms agree with Eq. (4.21) in Ref. [56]. Note that in that paper there are no poles 
in V s since the regularization used takes the Wilson lines off the light cone. 

In order to compute the anomalous dimension matrix for the reduced soft function, we 
should subtract the jet anomalous dimensions, according to Eq. (3.16), 



r /j (Pij, <**) = T/j (A • 0j,a a , e) - 5jj U k (wk, a a ,e) . 



(B.4) 



k=i 



The jet anomalous dimensions can be computed using Eq. (3.8) in the fundamental 
representation. At one loop this yields 



U k (w k ,a s ,e) = ^Cf 



In 



2(w fc • n k ) 



ni 



+ 0{a 2 s ). 



(B.5) 
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The subtracted terms are proportional to the unit matrix in colour space, so they affect only 
the diagonal elements of the anomalous dimension matrix. We end up with the following 
anomalous dimension for S, which is of course finite, 



a 



2 + \ In ( Pl2 P34) 



7T 



2 + \ In (p 13 p 24 ) 





+ 0(a 



It is straightforward to check that the general expressions in Eq. (5.16) and Eq. (5.6) indeed 
reduce to Eq. (B.6) and Eq. (B.3), respectively, upon evaluating the colour factors in the 
chosen basis and substituting the one-loop values for jk, <% and 5s- 
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